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Current interferometric gravitational wave detectors use test masses with mirror coatings formed 
from multiple layers of dielectric materials, most commonly alternating layers of SiCh (silica) and 
Ta2 05 (tantala). However, mechanical loss in the Ta2 0s/Si02 coatings may limit the design sensi- 
tivity for advanced detectors. We have investigated sources of mechanical loss in the Ta20s/Si02 
coatings, including loss associated with the coating-substrate interface, with the coating-layer in- 
terfaces, and with the bulk material. Our results indicate that the loss is associated with the bulk 
coating materials and that the loss of Ta2 0s is substantially larger than that of SiC>2. 

PACS numbers: 04.80.Nn, 95.55.Ym, 62.40.+i, 68.35.Gy 
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I. INTRODUCTION 

The sensitivity of designs for advanced interferomet- 
ric gravitational wave detectors, such as Advanced LIGO 
(Laser Interferometer Gravitational Wave Observatory), 
is limited in the frequency range from 10's to 100's of 
Hz by thermal noise from the main test masses and their 
suspensions pj. These test masses, which under current 
design plans will be formed from either fused silica or 
sapphire, will each be supported using a fused silica sus- 
pension, and will have multi-layer, dielectric mirror coat- 
ings |2|. In the current generation of gravitational wave 
detectors, the test mass mirror coatings are formed by 
ion-sputtering alternating layers of silicon dioxide (Si02) 
and tantalum pentoxide (Ta205). This type of coating 
was chosen because it can be made highly reflective in 
a narrow band around 1.064 /xm, the laser wavelength 
chosen for LIGO, while having very low absorption and 
scatter losses y|- 

For the past several years, many research groups in 
the gravitational wave field have been concerned that di- 
electric coatings could be an important source of ther- 
mal noise. This concern was strengthened by Levin's 
calculation |4| which indicated that mechanical losses 
in the mirror surface of a test mass could be much 
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greater than had been generally appreciated. Subse- 
quently, our investigations of mirror coatings showed that 
multi-layers of Ta205/Si02 and Ta20s/Al203 when ap- 
plied to fused silica substrates, add significant levels of 
mechanical loss 0, El- Using these results and models 
developed by Nakagawa Q, and Gretarsson @, we cal- 
culated that mechanical loss in the coatings would result 
in a level of thermal noise which would degrade the de- 
sign sensitivity of the planned advanced detectors by a 
significant amount. 

In the present study, we have performed a series 
of experiments to investigate the source of the loss in 
multi-layer tantala/silica coatings. We hypothesized that 
the loss would arise primarily either from the coating- 
substrate interface, from the coating-layer interfaces, 
from the bulk materials used in the coatings, or possi- 
bly from some combination of these sources. We have 
measured the loss of a series of coatings in which the 
number and thickness of the coating layers was chosen to 
test each of these dependencies. 



II. DESIGN OF COATING LOSS STUDY 
A. Experimental technique 

As in our previous experiments 5, 6], we determine the 
level of mechanical loss associated with a given coating 
by applying the coating to fused silica substrates and 
measuring the mechanical losses of a subset of resonant 
modes of the coated samples. 

The mechanical loss, <j), at a resonant frequency /o of 
a sample, is related to the quality factor, Q of the reso- 
nance by 4>(fo) = 1/Q. Assuming all other losses to be 
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FIG. 1: Coating on Thin sample with semicircular mask 
around weld. 

negligible, the total loss in a coated sample, </> C oated(/o), 
is equal to the sum of the intrinsic loss of the substrate 
plus any loss associated with the coating 0, |fj , 

Y^coatcd ^substrate V^coating (/o) (1) 

-^substrate 

where -Ecoating/^substrate is the ratio of energy stored in 
the coating to energy stored in the substrate. 

In this study, the substrates were 7.6 cm diameter fused 
silica disks of two different thicknesses, 2.5 cm and 0.25 
cm. The thinner disks had the advantage that the ef- 
fect from the coating was more pronounced; the coating 
loss had a greater contribution to the total loss because 
the energy ratio, boating /^substrate, was larger for the 
thin samples. However suspending these thin samples 
required the direct welding of a low loss silica suspen- 
sion to the edge of the fused silica substrates, see 0. 
To minimize any potential damage to the coatings by di- 
rect heat from welding, the coating was masked in a 1 
cm radius around the weld, see figure ^ I n addition, 
we used miniature welding torches that allowed the weld 
area to be sub-millimeter in scale. During welding, the 
glass Ri 1 cm from the weld remained only warm to the 
touch, which suggests that there was little to no impact 
on the coating. Although measurements of the loss fac- 
tors of the thicker samples were less sensitive to the ef- 
fects of coating losses, they had the advantage that they 
could be suspended using thin silk thread, see 

13, 

with- 
out any risk of physically altering the substrate or the 
coating. 

The use of both techniques allowed a more thorough 



investigation of mechanical losses in the coating over a 
wider range of frequencies @ . 

The substrates were made of Corning 7980 grade OA 
fused silica . Their faces were polished to sub-angstrom 
micro- roughness by Wave Precision Inc. ^(|> to emulate 
the required surface properties of actual test masses. Ex- 
cept where otherwise noted, all the coatings were applied 
by SMA/ Virgo 11] at llnstitut de Physique Nucleaire in 
Lyon France. 

To measure the quality factor of a given mode of a 
sample, we suspend the sample in vacuum and excite it 
to resonance using an electrostatic drive. We then remove 
the excitation signal and record the sample's motion as 
it freely rings down. The characteristic decay time, r, is 
the time required for the amplitude of motion to decrease 
by a factor l/e. The quality factor is then given by Q = 
ftfoT. The amplitude of the resonant motion was sensed 
interferometrically for the thick (2.5 cm) samples, and 
by polarimetry for the thin (0.25 cm) samples. At the 
vacuum pressure for the experiment, (P < 10~ 6 torr), 
air damping was negligible. Detailed descriptions of these 
ringdown techniques are given in [f| and [g respectively. 

We performed finite element analysis, with the pro- 
gram ALGOR [lj, to model the displacement of each 
mode of the samples. The relevant energy ratios for the 
coated samples were then calculated from the displace- 
ments 5] . If the intrinsic mechanical loss of the substrate, 
^substrate (/o)> is known (or is insignificant compared to 
the effective loss from the applied coating), the loss as- 
sociated with the coating may then be calculated using 
eqn. ^ 



Sample 


Total 


Optical Thickness 




Type 


Layers 


SiOa 


Ta 2 5 


Comments 


A 











Annealed only 


B 


2 


A/4 


A/4 


Coated & annealed 


C 


30 


A/4 


A/4 


Coated & annealed 


D 


60 


A/8 


A/8 


Coated & annealed 


E 


30 


A/8 


3A/8 


Coated & annealed 


F 


30 


3A/8 


A/8 


Coated & annealed 



TABLE I: Set of samples used to probe the source of the 
mechanical loss in the Ta20s/Si02 coating. 



B. Sequence of coatings studied 

We investigated three primary sources of mechanical 
loss in the coatings, which we postulated to be: 

• loss in the coating-substrate interface, 

• loss in the interfaces between the multiple coating 
layers, and 

• loss in the bulk coating materials. 




To investigate these hypotheses, we designed a set of 
coatings in which we varied the number and thickness 
of the coating layers to test for the three dependencies. 
Table[I]hsts the coatings and the substrate treatments in- 
vestigated. By comparing the measured mechanical loss 
for sample types B through F, we can test our three hy- 
potheses for the source of the coating loss. 

If the dominant source of coating mechanical loss is 
associated with the coating-substrate interface, then the 
total loss measured for sample types B-F would be ex- 
pected to be approximately equal. 

On the other hand, the loss may have originated pre- 
dominantly in the interfaces between the multiple coating 
layers. In that case the coating loss of sample type D, 
which has 60 layers, should be twice as large as the loss 
in sample types C, E and F which have 30 layers. 

Finally, the source of dissipation may be intrinsic to the 
bulk coating material. Under that scenario, the coating 
loss for sample types C, E, and F should vary with the 
abundances of the two coating materials. 

Clearly there also exists the possibility that the total 
mechanical loss of the coating is the sum of contributions 
from a combination of the above mechanisms. 

As part of the coating process, samples are heated to 
a few hundred degrees to reduce residual stresses in the 
coating. Previous experiments 

El III 111 have shown 
that heating uncoated fused silica samples can result in 
a significant decrease in loss. Coating Type A was not 
coated but was cleaned and heated to the same elevated 
temperature as the coated samples to ensure any change 
in substrate loss was taken into account in our analysis. 
We will refer to this process as "annealing" even though 
the samples were not raised to the standard annealing 
temperatures for fused silica. 



III. RESULTS 
A. Thin Sample Results 

We present below the results of our measurements on 
the thin samples. We acquired at least five measure- 
ments for each mode of a given sample configuration. As 
is common in Q experiments, we quote the result of the 
highest measurement. The reason for this selection pro- 
cess is that systematic errors, which in many experiments 
can equally increase or decrease one's measurement, will, 
with rare exception in our measurements, uniformly de- 
grade the Q. On the other hand, the statistical errors 
contribute a uniform, Gaussian background that intro- 
duces an error to the fit of Q. Data that does not exhibit 
uniform background noise is discarded. These statisti- 
cal errors are normally much smaller that the systematic 
fluctuations. Therefore we quote the largest value for our 
Q measurements with its associated statistical error. 

The suspension system, excitation system, and data 
analysis routines have all been tested to produce no lim- 
iting effect on measurements where the Q exceeds 80 mil- 




FIG. 2: Modes of the thin sample. The first mode (Butterfly) 
on the left and the second mode (Drumhead) on the right. 
The first row is face view and the second row is side view. 

lion. During these measurements, the dimensions of the 
suspension system and the position of the excitation sys- 
tem were altered with no observed change in the Q of 
the sample. Finally these measurements were performed 
at MIT and at Syracuse University with excellent agree- 
ment between the two facilities. 

The results are listed by mode in Table UTI The mode 
shapes are shown in Figure [3 The Butterfly- x and 
Butterfly-+ modes are the degenerate first mode with 
a typical resonant frequency of /i = 2.7 kHz. The addi- 
tional notation of "x" and "+" refer to the orientation of 
the nodal lines. The second mode, known as the drum- 
head mode, has a typical resonant frequency of fa = 4.1 
kHz. 

One immediate observation from the data is that 
for the coated samples, with the exception of the 2- 
layer coating, Qcoatcd << Quncoatod- Typical values for 
Quncoatcd indicate that the loss in the substrate is about 
2-3% of Qcoated and should not significantly affect the 
results for the coating loss. Moreover, we see that the 
annealing performed at the end of the coating run dra- 
matically increases the Q of an uncoated substrate. If 
we assume that the substrate in a coated sample also 
experiences the same increase in Q, then the loss in the 
substrate is reduced to about 1% of the loss in the coated 
sample. Thus, for the thin samples, we assume the loss 
measured in a coated sample is approximately the loss 
associated with the coating. 

, -^coating , /~n 

^coated ~ t-, ^coating \^ ) 

-^substrate 

To calculate the coating losses we require the fraction 
of the mode energy stored in the coating. This value is 
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Sample 


Butterfly x mode 


Butterfly + mode 


Drumhead mode 


Coating Thickness 


Type 


Number 


Quncoated 


Qcoated 


^coating 


Quncoated 


Qcoated 


^coating 


Quncoated 


Qcoated 


^coating 


Ssilica (M m ) 


Stantala (M m ) 


A 


1 


14.7 




N/A 


11.7 


43.6 


N/A 








Uncoated, Annealed at 600° C 


A 


2 


10.6 


42.4 


N/A 


13.9 


54.0 


N/A 








Uncoated, Annealed at 900° C 


B 


1 


24.2 


9.0 


2.4 


25.2 


8.0 


2.7 








0.183 


0.131 


B 


2 


6.32 


5.4 


4.0 


7.53 


6.5 


3.3 




6.4 


3.2 


0.183 


0.131 


C 


1 


18.4 


0.53 


2.7 


20.8 


0.55 


2.6 








2.75 


1.97 


C 


2 


23.1 


0.53 


2.7 


18.6 


0.54 


2.6 




0.43 


3.1 


2.75 


1.97 


D 


1 


17.3 


0.49 


2.9 


20.4 


0.55 


2.6 




0.44 


3.1 


2.75 


1.97 


D 


2 




0.54 


2.6 


43.6 


0.51 


2.8 








2.75 


1.97 


E 


1 


42.4 


0.40 


3.9 


54.0 


0.40 


3.9 




0.29 


5.1 


1.38 


2.95 


E 


2 


22.2 


0.40 


3.9 


20.3 


0.41 


3.8 








1.38 


2.95 


F 


1 




0.75 


1.8 




0.72 


1.8 








4.13 


0.983 


F 


2 




1.13 


1.2 




0.82 


1.6 




0.63 


2.0 


4.13 


0.983 



TABLE II: Results for thin sample measurements. All Q's are x 10 6 . All d>'s are x 10 



calculated from the displacement values obtained using 
finite element analysis package, ALG0R pjjj ■ This analysis 
yielded an energy ratio per coating thickness of 



dE/ds 

E 
dE/ds 

E 



1488 m" 



= 1575 m~ 



Butterfly mode 
Drumhead mode 



Then, given the coating thickness, s, we can rewrite 
the equation for (/> CO ating as 



^coated (/o) — s 



dE/ds 
E 



5 (/o) 



(3) 



The coatings are made such that the first layer applied 
to the substrate is a tantala layer. We convert a coat- 
ing's optical thickness to its physical thickness using the 
refractive indices n s iii ca = 1.45 and ntantaia = 2.03. 

Reviewing our results, listed in Table ITU we can test 
our hypotheses on the source of the coating loss. The 
first hypothesis was that the coating loss originated in the 
coating-substrate interface, in which case the loss would 
be independent of coating thickness. Were this idea true, 
then coating types B and C, which have 2 and 30 layers 
respectively, would have similar values for Qcoatod- In- 
stead, Qcoatod for these samples differ by a factor 10, 
and the coating loss, ^coating, which is scaled for coat- 
ing thickness, is roughly equal for the two coating types. 
Thus the loss is in the coating, not the coating-substrate 
interface. 

Our second idea was that the coating loss arose in the 
interfaces between the multiple coating layers. To test 
this hypothesis we can compare the loss for coating types 
C and D, which have the same total thickness but 30 
and 60 layers respectively. If the loss were dependent on 
the number of coating interfaces then <^ coa ting for coat- 
ing types C and D would differ by a factor 2. However, 




FIG. 3: The shapes of the four modes of the thick sample: 
Clover-4 or Butterfly (C4) in upper left; Fundamental Radial 
(F) in upper right; Asymmetric Drumhead (A) in lower left, 
and 2 nd Asymmetric (2A) in lower right. 



coating types C and D have similar values for </> C oating, 
indicating that the loss does not depend on the number 
of coating interfaces. 

Finally if the loss originates in the bulk coating mate- 
rial, then we would expect to see a dependence in ^ CO ating 
as the proportions of coating materials are varied in coat- 
ing types E, C, and F. Indeed boating does increase with 
increasing proportions of Ta20s, indicating that the loss 
is due to the bulk coating materials. 

To separate boating into the losses for the two coating 
materials, we must partition the energy in the coating 
into the amounts stored in its silica and tantala layers. 
Following an analysis by Landau and Lifshitz [Ti| , it can 
be shown that, for Poisson's ratio, a « 1, the energy 
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FIG. 4: The measured Q values for the four modes of each of the five thick samples with the coatings listed in Table Q] The 
color of the bar indicates whether the data was measured before (black) or after (white) the coating/ annealing process. 
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FIG. 5: A0(/o) plotted as a function of energy ratio for a 
typical sample. 
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FIG. 6: Same sample as in figure |H| except that A<^(/o) is 
corrected to account for the change in substrate loss due to 
annealing. 



stored in a coating, -E C oating <x Ys, where Y is Young's 
modulus and s is the coating thickness. For a multi-layer 
coating made up of two materials, we may write: 

A YlSl A Y2S2 A IA\ 

(/•coating = Tp <Pl + Tp <t>2 (4) 

-^coating ^coating ^coating ^coating 



where Yi, Si and 4>i are the Young's modulus, total thick- 
ness and loss angle of the i th coating material. We use 
the Young's moduli of 7.2 x 10 10 Pa for fused silica and 
1.4 x 10 11 Pa for tantalum pentoxide [r7j |. 

For a thin surface layer, the stress will be predomi- 
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nantly parallel to the surface. In this limit, the total 
Young's modulus for the multi-layer coating is given by 



•^coating -^coating 1 'I ~t~ 32*2 

Using eqn. \5\ we rewrite eqn.^Jas 

, siYi^i + s 2 Y 2 (f>2 



^coating 



siYi + s 2 Y 2 



(5) 



(6) 



Fitting the results for boating from coating types C, 
E, and F, we obtain the bulk loss values for silica and 
tantala 

Silica = (0.1 ± 0.1) x 1CT 4 (Butterfly mode) 
tantala = (4.7 ± 0.1) x 1(T 4 (Butterfly mode) 

^silica = (-0.1 ± 0.3) x 1CU 4 (Drumhead mode) 
^tantala = (6.1 ± 0.3) x 1CU 4 (Drumhead mode) 



As a final consistency check, we measured two sam- 
ples with coating type C that were coated by MLD [l8| . 



These measurements yielded 



^coating 



= 3.0 x 10" 4 for the 



butterfly mode, and </> CO ating = 3.2 x 1CU 4 for the drum- 
head and 2 nd butterfly modes. These results are in basic 
agreement with the samples coated at SMA/ Virgo 
and indicate that the coating loss is associated with the 
materials and not the coating process or manufacturer. 



B. Thick Sample Results 

1. Initial observations 

For the thick fused silica samples (7.6 cm x 2.5 cm 
thick), we measured the Q factors of 4 modes both be- 
fore and after each sample was coated or annealed. The 
results of these measurements are shown in figure^] The 
mode shapes, which were calculated using a finite element 
package [IJ , are shown in figure |31 

We suspended the thick samples in a loop of silk 
thread. Each sample was suspended multiple times, with 
the suspension length varied each time, and the highest 
Q factor for each mode was used in our analysis. Previ- 
ous experiments |0, 120I l2l) have shown that measured 
Q factors may be too low if the resonant frequency of a 
test mass happens to coincide with a resonant frequency 
of the suspension wires. Using this technique we have 
shown that suspension losses can be reduced to a negligi- 
ble level when measuring Q's of the order of a few times 
10 7 0. 

Before discussing the quantitative analysis of the re- 
sults, it is instructive to observe the trends in the data 
shown in figure 0] Firstly, it can be seen that the Q fac- 
tors for sample A, which was annealed, and for sample 
B, which had a 2-layer coating, were very similar, and 
that these Q factors were considerably higher than the 



other samples which had thicker coatings. This differ- 
ence indicates that the coating-substrate interface is not 
the dominant source of mechanical loss. 

Secondly, samples C and D, which have the same total 
coating thickness, but 30 and 60 coating layers respec- 
tively, have very similar Q factors. This result suggests 
that the individual multi-layer interfaces are not the dom- 
inant source of mechanical loss. 

Finally, the coatings on samples C, E, and F each 
have 30 layers and the same total thickness, but vary 
the proportions of SiC>2 and Ta20s. The measured Q's 
show that the loss increases with increasing proportions 
of Ta2C>5 in the coating. These results suggest that the 
bulk coating materials are the dominant source of me- 
chanical loss in the coatings, and that tantalum pentox- 
ide has a higher mechanical loss than silica. 

These observations are all consistent with the trends 
seen in the thinner samples discussed in section UlI Al 

It is important to note that the measured Q factors for 
the samples are mode-dependent. This dependence will 
be discussed in the following section. 



2. Quantitative analysis 

Calculating the coating loss for the thicker samples re- 
quires a slightly different method than for the thin sam- 
ples. For the thin samples, the loss in the substrate was 
negligible, so the coating loss could be calculated directly 
from the measured loss in the coated sample. For the 
thick samples, a much smaller fraction of the energy was 
stored in the coating. Thus the loss in the substrate must 
be considered when calculating the loss in the coating. 

For each mode of a thick coated sample, we calcu- 
late the energy ratio boating /^substrate- We then rewrite 
eqn. ^ as: 



A(p — '/'coated — ^substrate 



E, 



coating 



V^coating (7) 



To obtain ^coating) w e plot A<fi versus 
^coating/^substratc for each mode. Figure shows 
this plot for one sample with thirty alternating layers of 
Si02 and Ta20s. It can be seen that the data appears 
far from the expected straight line. 

From Figure 0] it can be seen that there were signif- 
icant changes in the Q factors of modes of a sample 
which was annealed but not coated. The Q's of some 
modes increased while others decreased. In our case, the 
mechanism for changes in Q is not known but we postu- 
late it may be due to some redistribution of stress in the 
samples, which may produce mode-dependent losses |l3j . 
Thus while using the measured values for the Q of the 
uncoated substrate seems invalid, replacing these values 
by the Q factors of the modes of the annealed, but un- 
coated mass, should in principle now allow us to obtain 
a value for C oating(/o)- The results of this are shown in 
figure El 
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It can be seen that there is a significant improvement 
of the fit to the data using this analysis. Carrying out 
the same types of analyses on the data from the other 
coated samples gave very similar results. 

3. Results 

The loss factors obtained for each type of coating stud- 
ied are shown in table Hill 



Coating 


Number of 


coating 


Type 


Samples 


B 


1 


(0.9 ±2.8) x 10" 


-4 


C 


2 


(2.7 ±0.7) x 10" 


-4 


D 


2 


(2.7 ±0.5) x 10" 


-4 


E 


2 


(3.7 ±0.5) x 10" 


-4 


F 


2 


(1.9 ±0.2) x 10" 


-4 



TABLE III: Mechanical loss factors for the various coating 
types calculated from the thick sample data. 

The results shown in tablc lllll confirm our earlier obser- 
vations indicating that the dominant source of mechani- 
cal loss is associated with the bulk coating materials, and 
also confirm that the tantalum pentoxide component of 
the coatings has a higher mechanical loss factor than the 
silica component. This last deduction can be treated 
more quantitatively. 

As with the thin samples, we use eqn.0and the results 
for ^coating from coating types C, E, and F, to calculate 
the loss for the silica and tantala layers. The results are: 

silica = (0.5 ± 0.3) x 1(T 4 

tantala = (4.4 ± 0.2) X 10~ 4 

We also performed a consistancy check on the thick 



samples by measuring two samples that had coatings 
of type C applied by MLD 18] . The measured coat- 
ing losses were C oating = 2.8 x 10~ 4 for sample 1, and 
'/'coating = 3.8 x 10~ 4 for sample 2. Like the thin sam- 
ples, the thick samples yielded similar results for the two 
coating manufacturers. In both cases the coatings from 
MLD showed slightly higher loss than the equivalent 
coatings produced by SM A/ Virgo [Til ]. 



IV. CONCLUSIONS 

For advanced interferometric gravitational wave 
observatories to be able to reach realistic astronomical 
distances, low mechanical loss mirror coatings are 
necessary. Observatories currently in operation use 
a multi-layer dielectric coating of silica and tantala. 
We have measured the loss in this coating to be 
4> = 2.7 x 10~ 4 , which is in basic agreement with the 
results from our previous work 0, Q . Those references 
discuss how this coating loss may effect the sensitivity of 
Advanced LIGO. We have shown in this work that the 
loss in this coating is predominantly the bulk loss in the 
tantala layers. Work is underway to investigate other 
coating materials and processes that may give lower 
mechanical loss while retaining the high reflectivity, low 
scatter, and low optical loss required by these advanced 
detectors. 



V. ADDENDUM 

The coating applied to the thick disk as reported in 
reference |(| is incorrect. The correct coating thickness 
of 4.66 (im results in a calculated coating </> of 5.2 x 10~ 4 . 



[1] LIGO Science Collaboration, P. Fritschel 
ed., Advanced LIGO System Design, 

LIGO-T010075-00-D (2001). Available at 
|http://www.ligo.caltech.edu/docs/T/T010075-00.pdf| 

[2] LIGO Science Collaboration, LSC White Paper on De- 
tector Research and Development (1999). Available at 
http://www.ligo.caltech.edU/docs/T/T990080-00.pdf 

[3] K. Srinivasan, D. Coyne, R. Vogt, LIGO docu- 
ment LIGO-T97016-00-D (1997) Available at 
|http://www.hgo.caltech.edu/docs/T/T970176-00.pdf| 

[4] Yu. Levin, Phys. Rev. D 57 659 (1998). 

[5] D. Crooks, P. Sneddon, G. Cagnoli, J. Hough, S. Rowan, 
M. M. Fejer, E. Gustafson, R. Route, N. Nakagawa, 
D. Coyne, G. M. Harry, A. M. Gretarsson, Class. Quan- 
tum Grav. 19 4229 (2002). 

[6] G. M. Harry, A. M. Gretarsson, P. R. Saulson, 
S. D. Penn, W. J. Startin, S. E. Kittelberger, 
D. R. M. Crooks, J. Hough, G. Cagnoli, N. Nakagawa, 



S. Rowan, M. M. Fejer, Classical and Quantum Gravity 
19 (2002) 897-918. 

[7] N. Nakagawa, A. M. Gretarsson, E. K. Gustafson, 
M. M. Fejer, Phys.Rev. D65 102001 (2002). 

[8] P. Sneddon and S. Rowan, private communication, pub- 
lication in preparation. 

[9] Corning Incorporated, Corning, New York 14831-0001, 
United States. 

[10] WavePrecisionlnc, 5390 Kazuko Court, Moorpark, CA 
93021, USA. 

[11] Service des Materiaux Avances/Virgo, Lyon, France. 
[12] ALGOR, Inc., 150 Beta Drive, Pittsburgh, PA 15238- 
2932 USA. 

[13] K. Numata, et al, Phys. Lett. A276 37 (2000). 

[14] B. S. Lunin and S. N. Torbin, Vestnik Moskovskogo Uni- 
versiteta, ser. 2, Khimia, 41 93 (2000). 

[15] S. D. Penn, G. M. Harry, A. M. Gretarsson, S. E. Kit- 
telberger, P. R. Saulson, J. J. Schiller, J. R. Smith, 



8 



S. O. Swords, Rev. Sci. Instrum. 72 3670 (2001). 

[16] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, 
3rd edition, Pergamon Press, Oxford, 1986. 

[17] P. J. Martin, A. Bendavid, et al., "Mechanical 
and Optical Properties of Thin Films of Tanta- 
lum Oxide Deposited by Ion-assisted Deposition", in 
"Thin films: Stresses and Mechanical Properties IV", 
Mat. Res. Soc. Symp. Proc. 308 583 (1993) 

[18] MLD Technologies, 2672 Bayshore Parkway, Suite #701, 



Mountain View, CA 94043 USA. 
[19] V. B. Braginksy, V. P. Mitrofanov, V. I. Panov, Systems 
with Small Dissipation, The University of Chicago Press, 
(1985) 

[20] J. E. Logan, N. A. Robertson and J. Hough, Phys. Lett. 

A170 352-358 (1992) 
[21] S. Rowan, S. M. Twyford, J. Hough, D. H. Gwo, 

R. Route, Phys. Lett. A246 471-8 (1998) 



